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We point out that the precision moasurernonts of the pscudo observablcs , Ab , and Apg performed 
at LEP and SLC suggest that in models with minimal-flavor-violation the sign of the Z-penguin 
amplitude is identical to the one present in the standard model. We determine the allowed range 
for the non-standard contribution to the Inami-Lim function C and show by analyzing possible 
scenarios with positive and negative interference of standard model and new physics contributions, 
that the derived bound holds in each given case. Finally, we derive lower and upper limits for the 
branching ratios of K'^ ■k'^vp, Kl ■k'^vv, Kl ji^ , B Xd.si'v, and Bd,s f-^ 1^^ 
within constrained minimal-flavor-violation making use of the wealth of available data collected at 
the Z-pole. 
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I. INTRODUCTION 

The effects of new heavy degrees of freedom appear- 
ing in extensions of the standard model (SM) can be 
accounted for at low energies in terms of effective oper- 
ators. The unprecedented accuracy reached by the elec- 
trowcak (EW) precision measurements performed at the 
high-energy e+e" colliders at LEP and SLC impose strin- 
gent constraints on the coefficients of the operators en- 
tering the EW sector. The best studied operators for 
constraining new physics (NP) are those arising from the 
vector boson two-point functions [1] , commonly referred 
to as oblique or universal corrections. A little less promi- 
nent are the specific left-handed (LH) contributions to 
the Zbb coupling [2], which are known as vertex or non- 
universal corrections. The tight experimental constraints 
[3] on the three universal parameters ei (T), 62 {U), and 
63 {S), and the single non-universal parameter (7b) 
pose serious challenges for any conceivable extension of 
the SM close to the EW scale. 

Other severe constraints concern extra sources of fla- 
vor and CP violation that represent a generic problem in 
many NP scenarios. In recent years great experimental 
progress has come primarily from the BaBar and Belle 
experiments running on the e+e" T(4S') resonance, 
leading not only to an impressive accuracy in the de- 
termination of the Cabibbo-Kobayashi-Maskawa (CKM) 
parameters [4] from the analysis of the unitarity trian- 
gle (UT) [5, 6], but also excluding the possibility of new 
generic flavor-violating couplings at the TeV scale. The 
most pessimistic yet experimentally well supported so- 
lution to the flavor puzzle is to assume that all flavor 
and CP violation is governed by the known structure of 
the SM Yukawa interactions. This assumption deflnes 
minimal-flavor-violation (MFV) [7-9] independently of 
the specific structure of the NP scenario [10]. In the case 
of a SM-like Higgs sector the resulting effective theory al- 
lows one to study correlations between K- and _B-dccays 
[10-12] since, by virtue of the large top quark Yukawa 



coupling, all flavor-changing effective operators involv- 
ing external down-type quarks are proportional to the 
same non-diagonal structure [10]. The absence of new 
CP phases in the quark sector does not bode well for a 
dynamical explanation of the observed baryon asymme- 
try of the universe. By extending the notion of MFV 
to the lepton sector [13], however, baryogenesis via lep- 
togenesis has been recently shown to provide a viable 
mechanism [14]. 

The purpose of this article is to point out that in 
MFV scenarios there exists a striking correlation be- 
tween the Z ^ bb pseudo observables (POs) R'l, Ab, and 
Ap'g measured at high-energy e+e~ colliders and all Z- 
pcnguin dominated low-energy flavor-changing-ncTitral- 
current (FCNC) processes, such as tt^vv, Kl 

TT°ui>, Kl -> /i+M", B Xd,svv, and B,ls l^'^fJ-' just 
to name a few.^ The crucial observation in this respect 
is that in MFV there is in general a intimate relation 
between the non-universal contributions to the anoma- 
lous Zbb couplings and the corrections to the flavor off- 
diagonal Zdjdi operators since, by construction, NP cou- 
ples dominantly to the third generation. In particular, 
all specific MFV models discussed in the following share 
the latter feature: the two-Higgs-doublet model (THDM) 
type I and II, the minimal-supcrsymmctric SM (MSSM) 
with MFV [8, 9], all for small tan/3, the minimal univer- 
sal extra dimension (mlJED) model [15], and the littlest 
Higgs model [16] with T-parity (LHT) [17] and degener- 
ate mirror fermions [18]. Note that we keep our focus on 
the LH contribution to the Z-penguin amplitudes, and 
thus restrict ourselves to the class of constrained MFV 
(CMFV) [11, 19] models, i.e., scenarios that involve no 
new effective operators besides those already present in 
the SM. As our general argument does not depend on the 



Of course, e'/^i Tr^i^"/ , B X^^sl^l and all exclusive 

b — > d{s)l^l~ transitions could be mentioned here too. 
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chirality of the new interactions it also applies to right- 
handed (RH) operators, though with the minor difficulty 
of the appearance of an additional universal parameter. 
Such an extension which covers large tan f3 contributions 
arising in a more general framework of MFV [10] is left 
for further study. 

This article is organized as follows. In the next sec- 
tion we give a model-independent argument based on the 
small momentum expansion of Feynman integrals that 
suggests that the differences between the values of the 
non-universal Zhb vertex form factors evaluated on-shell 
and at zero external momenta are small in NP models 
with extra heavy degrees of freedom. The results of the 
explicit calculations of the one-loop corrections to the 
non-universal LH contributions to the anomalous Zbb 
coupling in the CMFV models we examine confirm these 
considerations. They are presented in Sec. III. Sec. IV 
contains a numerical analysis of the allowed range for the 
non-standard contribution to the Z-penguin function C 
following from the presently available data. In this sec- 
tion also lower and upper bounds for the branching ratios 
of several rare K- and B-decays within CMFV based on 
these ranges are derived. Concluding remarks are given 
in Sec. V. Apps. A and B collects the analytic expressions 
for the non-universal contributions to the renormalized 
LH Zbb vertex functions in the considered CMFV mod- 
els and the numerical input parameters. 

II. GENERAL CONSIDERATIONS 

The possibility that new interactions unique to the 
third generation lead to a relation between the LH non- 
universal Zbb coupling and the LH flavor non-diagonal 
Zdjdi operators has been considered in a different con- 
text before [20] . Whereas the former structure is probed 
by the ratio of the width of the Z-boson decay into bot- 
tom quarks and the total hadronic width, the bot- 
tom quark left-right asymmetry parameter. At, and the 
forward-backward asymmetry for bottom quarks, Ap^, 
the latter ones appear in FCNC transitions involving Z- 
boson exchange. 

In the effective fleld theory framework of MFV [10] , one 
can easily see how the LH non-universal Zbb coupling and 
the LH flavor non-diagonal Zdjdi operators are linked 
together. The only relevant dimension-six contributions 
compatible with the flavor group of MFV stem from the 
SU{2) X U{1) invariant operators 

_ (1) 
04,2 - ^ {QlYuY^t'^j^Ql) (j)^T''D''<j), 

that are built out of the LH quark doublets Q l , the Higgs 
field (t>, the up-type Yukawa matrices Yu, and the SU{2) 
generators t". After EW symmetry breaking these op- 
erators are responsible both for the non-universal Zbb 
coupling (i = j = b) and the effective Zdjdi vertex 
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FIG. 1: Relative deviations (S„ for low (upper panel) and high 
values (lower panel) of M . The solid, dashed, and dotted curve 
correspond to n = 1,2, and 3, respectively. In obtaining the 
numerical values we have set Mz = 91 GeV and mt = 165 GeV. 
See text for details. 



{i 7^ i). Since all SM up- type quark Yukawa couplings 
y^i except the one of the top, yt, are small, one has 
{YiiY^)ji w yfVf*Vti so that only the top quark contri- 
bution to Eq. (1) matters in practice. 

That there exists a close relation is well-known in the 
case of the SM where the same Feynman diagrams re- 
sponsible for the enhanced top correction to the anoma- 
lous Zbb coupling also generate the Zdjdi operators. In 
fact, in the limit of infinite top quark mass the corre- 
sponding amplitudes are identical up to trivial CKM fac- 
tors. Yet there is a important difference between them. 
While for the physical Z ^ bb decay the diagrams are 
evaluated on-shell, in the case of the low-energy Z djdi 
transitions the amplitudes are Taylor-expanded up to ze- 
roth order in the off-shell external momenta before per- 
forming the loop integration. As far as the momentum 
of the Z-boson is concerned the two cases correspond to 
the distinct points = M| and = in phase-space. 

Observe that there is a notable difference between the 
small momentum expansion and the heavy top quark 
mass limit. In the former case one assumes <C -^wj i^t 
while in the latter case one has , ^ TOj . This dif- 
ference naturally affects the convergence behavior of the 
series expansions. While the heavy top quark mass ex- 
pansion converges slowly in the case of the non-universal 
one-loop SM corrections to the Zbb vertex [21], we will 
demonstrate that the small momentum expansion is well 
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behaved as long as the masses of the particles propagat- 
ing in the loop are not too small, i.e., in or above the 
hundred GeV range. 

The general features of the small momentum expansion 
of the one-loop Zbb vertex can be nicely illustrated with 
the following simple but educated example. Consider the 
scalar integral 

with p3 = 0. Note that we have set the space-time di- 
mension to four since the integral is finite and assumed 
without loss of generality ms 0. 

In the limit of vanishing bottom quark mass one has 
for the corresponding momenta Pi=P2 = 0- The small 
momentum expansion of the scalar integral Cq then takes 
the form 

oo / 2 \ ^ 

Co = Y.an(^) , (3) 

with = {pi — p-if = —2p\-pi. The expansion coeffi- 
cients a„ are given by [22] 

where 

1 /a;ilna;i a;2lna;2\ 

g(xi,X2) = , (5) 

a;i - X2 V 1 - 1 - a;2 / 

and Xi = raf/m^. Notice that in order to properly gen- 
erate the expansion coefficients a„ one has to keep Xi 
and X2 different even in the zero or equal mass case. The 
corresponding limits can only be taken at the end. 

In order to illustrate the convergence behavior of the 
small momentum expansion of the scalar integral in 
Eq. (3) for on-shell kinematics, we confine ourselves to 
the simplified case mi = m2 = M and ms = mt- We 
define 




for n = 1,2, The M-dependence of the relative de- 
viations 5n is displayed in Fig. 1. We see that while for 
values of M much below m* higher order terms in the 
small momentum expansion have to be included in or- 
der to approximate the exact on-shell result accurately, 
in the case of M larger than m* already the first correc- 
tion is small and higher order terms are negligible. For 
the two reference scales M = 80 GeV and M = 250 GeV 
one finds for the first three relative deviations 5^ numer- 
ically +9.3%, +1.4%, and +0.3%, and +1.1%, +0.02%, 
+0.00004%, respectively. 

It should be clear that the two reference points M = 
80 GeV and M = 250 GeV have been picked for a rea- 
son. While the former describes the situation in the SM, 



i.e.. the exchange of two pseudo Goldstone bosons and 
a top quark in the loop, the latter presents a possible 
NP contribution arising from diagrams containing two 
heavy scalar fields and a top quark. The above exam- 
ple indicates that the differences between the values of 
the non-universal Zbb vertex form factors evaluated on- 
shell and at zero external momenta are in general much 
less pronounced in models with extra heavy degrees of 
freedom than in the SM. In view of the fact that this 
difference amounts to a modest effect of around —30% in 
the SM [21] , it is suggestive to assume that the scaling of 
NP contributions to the non-universal parts of the Zbb 
vertex is in general below the ±10% level. This model- 
independent conclusion is well supported by the explicit 
calculations of the one-loop corrections to the specific 
LH contribution to the anomalous Zbb coupling in the 
CMFV versions of the THDM, the MSSM, the mUED, 
and the LHT model presented in the next section. 

We would like to stress that our general argument does 
not depend on the chirality of possible new interactions 
as it is solely based on the good convergence properties 
of the small momentum expansion of the relevant vertex 
form factors. Thus we expect it to hold in the case of RH 
operators as well. Notice that the assumption of MFV 
does not play any role in the flow of the argument itself 
as it is exerted only at the very end in order to establish 
a connection between the Zbb and Zdjdi vertices evalu- 
ated at zero external momenta by a proper replacement 
of CKM factors. Therefore it does not seem digressive 
to anticipate similar correlations between the flavor di- 
agonal and off-diagonal Z-penguin amplitudes in many 
beyond-MFV scenarios in which the modification of the 
flavor structure is known to be dominantly non-universal, 
i.e., connected to the third generation. See [23] for a se- 
lection of theoretically well-motivated realizations. These 
issues warrant a detailed study. 



III. MODEL CALCULATIONS 

The above considerations can be corroborated in an- 
other, yet model-dependent way by calculating explicitly 
the difference between the value of the LH Zdjdi vertex 
form factor evaluated on-shell and at zero external mo- 
menta. In the following this will be done in four of the 
most popular, consistent, and phenomenologically viable 
scenarios of CMFV, i.e., the THDM, the MSSM, both for 
small tan /3, the mUED, and the LHT model, the latter 
in the case of degenerate mirror fermions. All computa- 
tions have been performed in the on-shell scheme employ- 
ing the 't Hooft-Feynman gauge. The actual calculations 
were done with the help of the packages FeynArts [24] 
and FeynCalc [25], and LoopTools [26] and FF [27] for 
numerical evaluation. 
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FIG. 2: Examples of one-loop vertex diagrams that result in a 
non-universal correction to the Z djdi transition in assorted 
NP scenarios with CMFV. See text for details. 



Before presenting our results^ we collect a couple of 
definitions to set up our notation. In the limit of vanish- 
ing bottom quark mass, possible non-universal NP con- 
tributions to the renormalized LH off-shell Zdjdi vertex 
can be written as 



ji 



G 



_F _e_ 
/2 7r2 



Ml^V^)Vt^CM<l^)do,l^.d,,Z^ , (7) 



where i = j = b and i ^ j in the flavor diagonal and 
off-diagonal cases. Gp, e, and Cw denote the Fermi 
constant, the electromagnetic coupling constant, the sine 
and cosine of the weak mixing angle, respectively, while 
Vij are the corresponding CKM matrix elements and the 
subscript L indicates that the interactions involve LH 
down-type quark fields only. 

As a measure of the relative difference between the 
complex valued form factor CNp(g^) evaluated on- shell 
and at zero momentum we introduce 



(5Cnp = 1 — 



ReCNp(g" =0) 
RcCNp(g2 = Ml) 



(8) 



In the THDM with vanishing tree-level FCNCs, the 
only additional contribution to the Z djdi transi- 
tions with respect to the SM comes from loops containing 
charged Higgs bosons, H^, and top quarks, t. An exam- 
ple of such a contribution is shown on the top left-hand 
side of Fig. 2. The correction depends on the mass of 
the charged Higgs boson, M^, and on the ratio of the 
vacuum expectation value of the Higgs doublets, tan (3. 
Models of type I and II differ in the way quarks couple 
to the Higgs doublets: in the type I scenario both the 
masses of down- and up-type quarks are generated by 
one of the doublets, like in the SM, while in the type II 



theory one of the doublets generates the down-type and 
the second one generates the up-type masses, like in the 
MSSM. In our case only the coupling to the top quark is 
relevant, so that we do not need to actually distinguish 
between types I and II. 

To find (5Cthdm we have computed analytically the 
one-loop charged Higgs corrections to Eq. (7) reproducing 
the result of [28]. The analytic expression for CTHDM('i''^) 
can be found in Eq. (Al). The dependence of i^Ctdhm 
on can be seen in the first panel of Fig. 3. The red 
(gray) band underlying the solid black curve shows the 
part of the parameter space satisfying the lower bound 
M| > 295 GeV following from B X^-f in the THDM 
of type II using the most recent SM prediction [29] . This 
tan /3 independent bound is much stronger than the one 
from the direct searches at LEP corresponding to > 
78.6 GeV [30], and than the indirect lower limits from a 
number of other processes. In model I, the most impor- 
tant constraint on comes from i?" [31]. As the cor- 
responding bound depends strongly on tan/3 we do not 
include it in the plot. While the decoupling of SGthdm 
occurs slowly, we find that the maximal allowed rela- 
tive suppression of ReCTHDM(9^ = M'^) with respect to 
Re Cthdm('Z^ = 0) is below 2% and independent of tan /3, 
as the latter dependence exactly cancels out in Eq. (8) . 
In obtaining the numerical values for 5Ctdhm we have 
employed M„ ^ 80 GeV, = 91 GeV, rrit = 165 GeV, 
and = 0.23. If not stated otherwise, the same numeri- 
cal values will be used in the remainder of this article. We 
assess the smallness of (5Cthdm as a first clear evidence 
for the correctness of our general considerations. 

In the case of the MSSM with conserved i?-parity, we 
focus on the most general realization of MFV compatible 
with renormalization group (RG) invariance [10]. In this 
scenario CKM-type flavor- and CP-violating terms ap- 
pear necessarily in the down- and up-type squark mass- 
squared matrices due to the symmetry principle underly- 
ing the MFV hypothesis. The explicit form of the phys- 
ical up-type squark mass matrix used in our analysis is 
given in Eq. (A8). We assume universality of soft super- 
symmetry (SUSY) breaking masses and proportionality 
of trilinear terms at the EW scale, ^ so that neutralino and 
gluino contributions to flavor-changing Z —^ djdi transi- 
tions are absent. This additional assumption about the 
structure of the soft breaking terms in the squark sector 
has a negligible effect on the considered FCNC processes 
[32]. Moreover, in the small tan/3 regime both neu- 
tralino and neutral Higgs corrections to Z ^ bb turn out 



The analytic expressions for the renormaUzed Zdj di vertex func- 
tions in the considered CMFV models are collected in App. A. 



If universality of soft SUSY breaking masses and proportional- 
ity of trilinear terms is assumed at some high-energy scale ofE- 
diagonal entries are generated by the RG running down to the 
EW scale. We ignore this possibility here. 

In [33] it has been pointed out that in scenarios characterized by 
large values of the higgsino mass parameter, i.e., 1 TeV, 

the MFV MSSM with small tan /3 is not necessarily CMFV due 
to the presence of non-negligible gluino corrections in AB = 2 
amplitudes. This observation is irrelevant for our further discus- 
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to be insignificant [34]. Therefore only SUSY diagrams 
involving chargino, , and stop, t, exchange are relevant 
here. An example of such a contribution can be seen on 
the top right side of Fig. 2. A noticeable feature in the 
chosen setting is that large left-right mixing can occur in 
the stop sector, leading to both a relatively heavy Higgs 
in the range 120 GeV < < 135 GeV and a stop mass 
eigenstate, say, ti, possibly much lighter than the remain- 
ing squarks. Such a scenario corresponds to the "golden 
region" of the MSSM, where all experimental constraints 
are satisfied and fine-tuning is minimized [35]. For what 
concerns the other sfermions we neglect left-right mixing 
and assume that all squarks and sleptons have a common 
mass rriq and nii, respectively.^ 

In order to find the complete MSSM correction 
(^Cmssm we have calculated analytically the one-loop 
chargino-up-squark corrections to Eq. (7) and combined 
it with the charged Higgs contribution. Our result for 
C^±{q^) agrees with the one of [34]^ and is given in 
Eq. (A2) . The region of parameters in which the SUSY 
corrections to the LH Zdjdi vertices are maximal cor- 
responds to the case of a light stop and chargino. In 
our numerical analysis we therefore focus on these sce- 
narios. We allow the relevant MSSM parameters to float 
freely in the ranges 2 < tan/3 < 6, |^| < 500 GeV, and 
M < ITeV for M = M^, M2,mi^,mq,mi. The value 
of the trilinear coupling A„ is computed from each ran- 
domly chosen set of parameters /i, mf_^, and rriq. The 
calculation of 6 ^ 57 and b sl'^l' that is used to con- 
strain the parameter space introduces also a dependence 
on the gluino mass Mg. We choose to vary Mg in the 
range 241 GeV < Mg < ITeV [36]. 

The MSSM parameter space is subject to severe ex- 
perimental and theoretical constraints. We take into ac- 
count the following lower bounds on the particle masses 
[30]: M± > 78.6 GeV, M|; > 94 GeV, m^^ > 95.7 GeV, 
rriq > 99.5 GeV, and > 73 GeV. In the considered 
parameter space the requirement of the absence of color 
and/or charge breaking minima sets a strong upper limit 
of around 3TeV on the absolute value of [37]. As 
far as the lightest neutral Higgs boson is concerned, we 
ensure that M° > 114.4 GeV [3], including the dominant 
radiative corrections [38] to its tree-level mass. Further 
restrictions that we impose on the SUSY parameter space 
are the p parameter [39], R° [34], and the inclusive b ^ sj 
[40] and b sl~^l~ [41] branching fractions. To find the 
boundaries of the allowed parameter space we perform 
an adaptive scan of the eight SUSY variables employing 
the method advocated in [42]. 

The dependence of (5Cmssm on the lighter chargino 



^ A strict equality of left-handed squark masses is not allowed due 
to the different D-terms in the down- and up-type squark sector. 



For our purposes this difference is immaterial. 

The last equation in this article has a typographic error. 

Passarino-Veltman function Cn should read C12. 
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mass M^^ is illustrated in the second plot of Fig. 3. Re- 
gions in the M^_^~ SCmssm plane where the absolute value 
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of the correction Re Cmssm (9^ = 0) amounts to at least 
2%, 4%, 6%, and 10% of the SM value ReCsM(g^ = 0) 
are indicated by the red (gray), green (light gray), blue 
(dark gray), and black points, respectively. No con- 
straints are imposed for the black points while the col- 
ored (grayish) ones pass all the collider and low-energy 
constraints mentioned above. 

Three features of our numerical explorations deserve 
special mention. First, the maximal allowed relative 
size of the correction ReCMssM('?^ — 0) amounts to 
less than of ReCsM('?^ = 0). Second, the magni- 
tude of the possible deviation (5 Cmssm is strongly anti- 
correlated with the absolute size of Re Cmssm (q^ = 0). 
While small corrections Re Cmssm (9^ = 0) allow for large 
values of (5Cmssm the latter difference decreases rapidly 
with increasing Re Cmssm (?^ = 0). Third, the correction 
'^C'mssm decouples quickly for heavy charginos. These 
features imply that ^Cmssm is small if one requires (i) 
the relative size of ReCMssM(9^ = 0) to be observable, 
i.e., to be bigger than the SM uncertainty of the univer- 
sal Z-penguin function^ and (ii) the chargino mass M^^ 
to be not too light. For example, all allowed points sat- 
isfy |(5Cmssm| < 0.01 if one demands |ReCMSSM('7^ = 
0)/ReCsM(g^ = 0)1 > 0.05 and M|; > 300 GeV. On the 
other hand, if the masses of the lighter chargino and stop 
both lie in the hundred GeV range, (5 Cmssm frequently 
turns out to be larger than one would expected on the 
basis of our model-independent considerations. 

The large corrections (5Cmssm can be traced back to 
the peculiar structure of the form factor C^± {q"^). While 
in the limit of vanishing external .^-boson momentum 
the first three terms in Eq. (A2) all approach a con- 
stant value the fourth one scales like ^^/-^Iusy with 
■^SUSY = iiiiii(Af , rrij^ ). Naively, one thus would ex- 
pect the general argument given in the last section to 
hold. Yet for large left-right mixing in the stop sec- 
tor, which permits a relatively heavy Higgs mass of 
> 120 GeV, it turns out that the first three con- 
tributions tend to cancel each other and, in turn, the 
size of (5Cmssm is controlled by the fourth term. Then 
•^Cmssm oc M'^/M^-^jgy and the correction i5Cmssm can 
be sizable if Msusy is close to the EW scale. The ob- 
served numerical cancellation also explains why (5Cmssm 
is typically large if Re Cmssm (s^ = 0) is small and vice 
versa. It should be clear, however, that the large devia- 
tion 5Cmssm are ultimately no cause of concern, because 
|ReCMSSM(9^ = 0)/ReCsM(<?^ = 0)| itself is always be- 
low 10%. In consequence, the model-independent bound 
on the NP contribution to the universal Z-penguin func- 
tion that we will derive in the next section does hold in 
the case of the CMFV MSSM. 



Among the most popular non-SUSY models in ques- 
tion is the model of Appelquist, Cheng, and Dobrescu 
(ACD) [15]. In the ACD framework the SM is extended 
from four-dimensional Minkowski space-time to five di- 
mensions and the extra space dimension is compactified 
on the orbifold jZi in order to obtain chiral fermions 
in four dimensions. The five-dimensional fields can equiv- 
alently be described in a four-dimensional Lagrangian 
with heavy Kaluza-Klein (KK) states for every field that 
lives in the fifth dimension or bulk. In the ACD model 
all SM fields are promoted to the bulk and in order to 
avoid large FCNCs tree-level boundary fields and inter- 
actions are assumed to vanish at the cut-off scale. ^ A 
remnant of the translational symmetry after compactifi- 
cation leads to KK-parity. This property implies, that 
KK states can only be pair-produced, that their virtual 
effect comes only from loops, and causes the lightest KK 
particle to be stable, therefore providing a viable dark 
matter candidate [46]. 

In the following we will assume vanishing boundary 

terms at the cut-off scale and that the ultraviolet (UV) 
completion does not introduce additional sources of fia- 
vor and CP violation beyond the ones already present 
in the model. These additional assumptions define the 
mUED model which then belongs to the class of CMFV 
scenarios. The one-loop correction to F™^^^ is found 
from diagrams containing apart from the ordinary SM 
fields, infinite towers of the KK modes corresponding to 
the T4^-boson, VT^^^, the pseudo Goldstone boson, 

the SU{2) quark doublets, Qg(fc), and the SU{2) quark 
singlets, Uq(^k)- Additionally, there appears a charged 
scalar, a^^, which has no counterpart in the SM. A pos- 
sible diagram involving such a KK excitation is shown on 
the lower left side in Fig. 2. Since at leading order the 
Z djdi amplitude turns out to be cut-off independent 
the only additional parameter entering F^^^^ relative to 
the SM is the inverse of the compactification radius 1/R. 
The analytic expression for CmUED(9^) can be found in 
Eq. (A9). 

For a hght Higgs mass of = 115 GeV a careful 
analysis of oblique corrections [47] gives a lower bound 
of 1/i? > 600 GeV, well above current collider limits [48]. 
With increasing Higgs mass this constraint relaxes sig- 
nificantly leading to 1/i? > 300 GeV [47, 49]. Other 
constraints on 1/i? that derive from i?^ [50], the muon 
anomalous magnetic moment [51], and flavor observables 
[52-54] are in general weaker. An exception is the inclu- 
sive B Xg^/ branching ratio. Since the SM prediction 
[29] is now lower than the experimental world average by 
more than la and the one-loop KK contributions inter- 
fere destructively with the SM 6 — ^ 57 amplitude [53, 55], 



The overall uncertainty of ReCsM(9^ = 0) amounts to around 
±3%. It is in equal shares due to the parametric error on the 
top quark mass, the matching scale uncertainty in the next-to- 
leading order result [43], and two-loop EW effects that are only 
partly known [44]. 



^ Boundary terms arise radiatively [45] . They effect the Z ^ dj di 
amplitude first at the two-loop level. Since we perform a leading 
order analysis in the ACD model its consistent to neglect these 
effects. 
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B Xs"f provides at leading order the lower bound 
1 /R > 600 GeV independent from the Higgs mass [56] . 
The 1 /R dependence of (5CmUED is displayed in the third 
plot of Fig. 3. In the range of allowed compactification 
scales, indicated by the red (gray) stripe, the suppression 
of ReCmUEDC?^ ~ M^) compared to ReCi„uED(9^ = 0) 
amounts to less than 5%, the exact value being almost 
independent of 1/R. This lends further support to the 
conclusion drawn in the last section. We finally note that 
our new result for Ci„ued(9^) coincides for = with 
the one-loop KK contribution to the Z-penguin function 
calculated in [52]. 

Another phenomenologically very promising NP sce- 
nario is the LHT model. Here the Higgs is a pseudo Gold- 
stone boson arising from the spontaneous breaking of an 
approximate global SU (5) symmetry down to SO{5) [16] 
at a scale /. To make the existence of new particle in 
the 1 TeV range consistent with precision EW data, an 
additional discrete Z2 symmetry called T-parity [17], is 
introduced, which as one characteristic forbids tree-level 
couplings that violate custodial SU (2) symmetry. In the 
fermionic sector, bounds on four fermion operators de- 
mand for a consistent implementation of this reflection 
symmetry the existence of a copy of all SM fermions, 
aptly dubbed mirror fermions [18]. The theoretical con- 
cept of T-parity and its experimental implications resem- 
ble the one of i?-parity in SUSY and KK-parity in uni- 
versal extra dimensional theories. 

Unless their masses are exactly degenerate, the pres- 
ence of mirror quarks leads in general to new flavor- and 
CP-violating interactions. In order to maintain CMFV 

we are thus forced to assume such a degeneracy here. In 
this case contributions from particles that are odd under 
T-parity vanish due to the Glashow-Iliopoulos-Maiani 
(GIM) mechanism [57], and the only new particle that 
affects the Z djdi transition in a non-universal way 
is a T-even heavy top, T+. A sample diagram involving 
such a heavy top, its also T-even partner, i.e., the top 
quark t, and a pseudo Goldstone field, G^, is shown on 
the lower right-hand side of Fig. 2. In turn, F^/^^ de- 
pends only on the mass of the heavy quark T+, which 
is controlled by the size of the top Yukawa coupling, by 
/, and the dimensionless parameter xl = Aj/(Af + A|). 
Here Ai is the Yukawa coupling between t and T+ and A2 
parametrizes the mass term of T+. In the fourth panel of 
Fig. 3, we show from bottom to top (5Clht as a function 
of XL for / = 1,1-5, and 2 TeV. The colored (grayish) 
bands underlying the solid black curves correspond to 
the allowed regions in parameter space after applying the 
constraints following from precision EW data [58] . As NP 
effects in the quark flavor sector of the LHT model with 
CMFV are generically small [59, 60], they essentially do 
not lead to any restrictions. We find that the maximal 
allowed suppression of ReCLHT('i'^ = ^^1) with respect 
to ReCLHT(g^ = 0) is slightly bigger than 3%. This 
feature again confirms our general considerations. Our 
new result for Clht('Z^) given in Eq. (AlO) resembles for 
= the analytic expression of the one-loop correction 



to the low-energy Z-penguin function calculated in [60]. 
Taking into account that the latter result corresponds to 
unitary gauge while we work in 't Hooft-Feynman gauge 
is essential for this comparison. In particular, in our case 
no UV divergences remain after GIM, as expected on 
general grounds [61]. 

At this point a further comment concerning gauge in- 
variance is in order. It is well known that only a proper 
arrangement of, say, e+e^ ff, including all contri- 
butions related to the Z-boson, purely EW boxes, and 
the photon, is gauge invariant at a given order in pertur- 
bation theory. In flavor physics such a gauge indepen- 
dent decomposition [62] is provided by the combinations 
X = C + B"'', Y = C + B^*^~, and Z = C + D/4 
of Inami-Lim functions [63]. Given the normalization 
of Eq. (7), NP contributions to the universal Z-penguin 
function C are characterized by ReC^p{q'^ = 0) in our 
notation, while B"^" and i?' ' represent the contribu- 
tion of EW boxes with neutrino and charged lepton pairs 
in the flnal state. D stems from the off-shell part of 
the magnetic photon penguin amplitude. Since we want 
to relate in a model-independent way observables de- 
rived from e+e~ — > // to observables connected with 
the di djvv and di djl^l~ transitions, we also have 
to worry about the potential size of corrections that are 
not associated with the Z-boson. 

At the Z-pole, the total cross-section of e+e" ff 
is completely dominated by Z-boson exchange. While 
purely EW boxes are vanishingly small, the bulk of the 
radiative corrections necessary to interpret the measure- 
ments are QED effects. It is important to realize that 
these QED corrections are essentially independent of 
the EW ones, and therefore allow the anomalous Zbb 
couplings to be extracted from the data in a model- 
independent manner. Certain SM assumptions are nev- 
ertheless employed when extracting and interpreting the 
couplings, but considerable effort [3] has been expended 
to make the extraction of the POs i?", Ab, and Apg as 
model-independent as possible, so that the meanings of 
theory and experiment remain distinct. 

In the case of the di djuv and — » djl~^l~ ob- 
servables theoretical assumptions about the size of the 
EW boxes are unfortunately indispensable. Our explicit 

analysis of the considered CMFV models reveals the fol- 
lowing picture. In the THDM, the NP contributions 
AB"^ = B"^-Bll^ and AS'"^'" = S'"^'" -B^^f vanish 
identical [64], while their relative sizes compared to the 
corresponding SM contributions amount to at most 
and ^5*% in the MSSM [65] and less than +1% in both 
the mUED scenario [52] and the LHT model [60]. The 
numbers for AB^^ , AB'^' , and AC quoted here and 
in the following refer to the 't Hooft-Feynman gauge. 
Moreover, contributions to the EW boxes are found to 
be generically suppressed by at least two inverse powers 
of the scale of NP using naive dimensional analysis [66]. 
In view of this, the possibility of substantial CMFV con- 
tributions to the EW boxes seems rather unlikely. The 
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actual size of the NP contribution AD = D — Dsm to 
the ofF-shell magnetic photon penguin function D has es- 
sentially no impact on our conclusions. The treatment of 
AB"", AS' and AD in our numerical analysis will 
be discussed in the next section. 



IV. NUMERICAL ANALYSIS 

Our numerical analysis consists of three steps. First we 
determine the CKM parameters p and f/ from an analysis 
of the universal UT [9] The actual analysis is performed 
with a customized version of the CKMfitter package [5]. 
Using the numerical values of the experimental and the- 
oretical parameters collected in App. B we find 

p = 0.160 ±0.031, 77 = 0.326 ± 0.012 . (9) 

The given central values are highly independent of 

rrit, but depend mildly on the hadronic parameter 

^ = {fB,B]^^)/{fBaB]^^) determined in lattice QCD. 
Since in our approach theoretical parameter ranges are 
scanned, the quoted 68% confidence levels (CLs) should 
be understood as lower bounds, i.e., the range in which 
the quantity in question lies with a probability of at least 
68%. The same applies to all CLs and probability regions 
given subsequently. 

In the second step, we determine the allowed ranges 
of AC and the NP contribution ACf = Cf{mb) - 
Cj%-^.i{inb) to the effective on-shell magnetic photon pen- 
guin function from a careful combination of the results of 
the POs i?^, Ab, and Ap^ [3] with the measurements of 
the branching ratios oi B ^ Xg-y [68] and B Xsl~^l~ 
[69]. In contrast to [12], we do not include the available 
experimental information on — » ir'^uD [70] in our 
global fit, as the constraining power of the latter mea- 
surement depends in a non-negligible way on how the 
experimental CL of the signal [71] is implemented in the 
analysis. 

Third, and finally, we use the derived ranges for the 
Inami-Lim functions in question to find lower and upper 
bounds for the branching ratios of the rare decays 
n^vD, Kl TT^vi/, Kl M"*^/^") B X^^sUU, and 
Bd,s Ai^Ai" within CMFV. 

Our data set includes all POs measured at LEP and 
SLC that are related to the Z ^ bb decay. It is given in 
Tab. I. Concerning the used data we recall that the ratio 
4/3 Apg/^e is lower than the direct measurement of Ah 
by 1.6(7, and lower than the SM expectation for Ab by 
3.2 cr [3]. Whether this is an experimental problem, an 
extreme statistical fluctuation or a real NP effect in the 



^ If the unitarity of the 3x3 CKM matrix is relaxed sizable devi- 
ations from Vth ~ 1 arc possible [67]. Wc will not consider this 
possibility here since it is not covered by the MFV hypothesis 
which requires that all flavor and CP violation is determined by 
the structure of the ordinary 3x3 SM Yukawa couplings [10]. 



TABLE I: Results and correlations for the Z -» 66 POs of the 
fit to the LEP and SLC heavy flavor data taken from [3]. 



Observable 


Result 


Rb 


Ab 


40,6 

^FB 


K 

Ab 
40,6 

^FB 


0.21629 ± 0.00066 
0.923 ±0.020 
0.0992 ± 0.0016 


1.00 


-0.08 
1.00 


-0.10 
0.06 
1.00 



bottom quark couplings is up to date unresolved. " In 

fact, the relative experimental error in is much larger 
than the ones in the total Z ^ bb rate, and Ab, where 
no anomalies arc observed. Furthermore, the extracted 
value of the anomalous LH coupling of the bottom quark 
agrees with its SM value because of the strong constraint 
given by i?". This strong constraint carries over to our 
results, which do not depend notably on whether Ap^ 
is included in or excluded from the data set. We as- 
sume that statistical fluctuations are responsible for the 
observed discrepancy and include A^^ in our global fit. 

The actual calculations of B!^, Ab, and Apg used in our 
analysis are performed with ZFITTER [75] , which includes 
the SM purely EW, QED and QCD radiative effects, pho- 
ton exchange and 7-Z interference that are necessary to 
extract the POs in a model-independent manner. For 
the purpose of our analysis, ZFITTER has been modified 
to include possible NP contributions to the Zbb vertex in 
the parametrization of Eq. (7). The Higgs mass is allowed 
to vary freely in the range 100 GeV < < 600 GeV. 
Since i?° is largely insensitive to the mass of the Higgs 
boson this conservative range has no noticeable impact 
on our results. The input values of the other parameters 
entering Ab, and Ap^ are collected in App. B. 

The experimental results that we consider in connec- 
tion with B Xg^ and B Xsl~^l~ are summarized 
in Tab. II. The given weighted average of the branc;h- 
ing ratio B{B — !■ Xg-y) corresponds to a photon energy 
cut Ej > 1.6 GeV in the B-meson rest frame, while for 
B{B Xsl^l~) the experimental data in the low-g^ 
region 1 GeV^ < (7^ < 6 GeV^ of the dilepton invari- 
ant mass squared, averaged over electrons and muons 
are shown. Our calculations rely on [29] in the case 
of B — > and on [77] in the case of B — > Xsl~^l~. 



It has been known for some time that measured with respect 
to thrust axis is not infrared (IR) safe [72]. Recently, an IR 
safe definition of has been suggested [73] which defines the 
direction of the asymmetry by the jet axis after clustering the 
event with an IR safe flavor jet-algorithm [74]. Given the long- 
standing discrepancy in it would be interesting to reanalyze 
the existing data using this alternative definition. 
The default flags of ZFITTER version 6.42 are used, except for set- 
ting ALEM = 2 to take into account the externally supplied value 
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1-CL 



TABLE II: World averages of I3{B X^j) for > 1.6 GeV 

and B{B Xsl+T) for 1 GeV^ < < 6GeV^ 



Observable 


Result 


B{B X,7) X 10* 


3.55 ±0.26 [76] 


B{B XJ+r ) X 10'' 


1.60 ±0.51 [69] 



The used numerical input parameters can be found in 
App.B. Unlike [12], we do not include B Xgl^l^ 
data on the regions 0.04 GeV^ < < 1 GeV^ and 
14.4 GeV^ < < 25GeV^ in our analysis. The rea- 
son for this omission is twofold. First, in these regions 
the differential B Xsl~^l~ rate is less sensitive to AC 
than in the low-q^ region. Second, for high the the- 
oretical uncertainties are larger with respect to the ones 
that affect the low-g^ region. An inclusion of the lat- 
ter two constraints would therefore make the fit more 
complicated, but it would not improve the quality of the 
obtained results. 

Before we present our final results, additional com- 
ments on the used methodology concerning AB'^'^, 



AB' ' , AC, AD, and AC f are in order. We begin with 
ACf which enters both B ^ Xs^f and B XJ+l^ . A 
well-known way to avoid the B Xgj constraint con- 
sists in having a large positive NP contribution AC^^ 
that approximately reverses the sign of the amplitude 
A{b 57) cx C^^{mb) with respect to the SM and leaves 
B{B Xs-f) oc \Cf{mb)\^ unaltered within experimen- 
tal and theoretical uncertainties. In our analysis, we add 
ACy^ to the top quark contribution of the SM, keeping 
TOf, that multiplies this combination renormalized at m^. 
This rescaling is motivated by the observation [78], that 
in this way most of the logarithmic enhanced QCD cor- 
rections are properly taken into account. We recall that 
C^^-^{mb) ~ —0.38 in this approach. 

Both the value and the sign of Cj^{mt,) play an impor- 
tant role in the B — + Xsl^l~ decay rate [79]. By contrast 
the dependence of B{B — > X^Z+Z") on D is relatively 
weak. Nevertheless, for suitable chosen values of AD 
the B Xsl^l~ constraint can be always satisfied even 
in the case of the non-SM solution of ACf^. In conse- 
quence, AD is not well constrained by the data used, and 
we decided to scan AZ? in the range ±1 for the best fit 
value. This choice is rather generous since in the CMFV 
scenarios that we consider one has \AD\ < JDsmI with 
DsM — —0.49 throughout the allowed parameter space 
[52, 60, 65, 80]. We verified that even larger variations 
have basically no effect on the extraction of the allowed 
range for AC, since the Z ^ bb POs do not depend on 
AD. The impact of AD on the bounds of ACf will be 
discussed below. 

Precision data on Ab, and Ap^ lead to a tight, 
highly model-independent constraint on Re Cnp (9^ = 
Af^). The allowed range of AC can then be calculated 
from the identity AC = (I + 5CNp)Re CNp(g^ = M|) in 
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FIG. 4: The upper (lower) panel displays the constraints on 

ACf and AC within CMFV scanning AB' ' in the range 
±0.1 (set to zero) that follow from a combination of the Z ^ bb 
POs with flavor observables. The colors encode the frequentist 
1 — CL level and the corresponding 68% and 95% probability 
regions as indicated by the bars on the right side of the panels. 
See text for details. 



any given model of NP where SC^p is known. To carry 
out the analysis in a generic way, one, however, needs 
to make an assumption about the size of 5Cnp . Guided 
by the results of Sees. II and III we allow JCnp to float 
in the range ±0.1. We note that larger variations with, 
say, an absolute value of [iJCsmI — 0.3, still lead to the 
conclusion that large negative values of AC that would 
reverse the sign of Csm — 0.78 are highly disfavored. 

The only EW box that enters the determination of 
ACf and AC in our case is AB^ ' . To explore the 
impact of the size of EW boxes on the fit results we 
consider two scenarios. In the first we allow Ai?' ' to 
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vary in the range ±0.1 while in the second we assume 



AS 



0. The former choice seems conservative, as 



relative to the SM value ~ 0.18 possible AS' 

contributions amount to only ^5^% in the MSSM with 
MFV and small tan /3 [65] and to below +1% in both the 
mUED [52] and the LHT model with degenerate mir- 
ror fermions [60]. In fact, the actual size of AB' ' , 
which enters our fit through B{B — » XJ~^l~), does only 
have a marginal effect on the results, because C is al- 
ready tightly constraint by the combination of At, 
and Ap'^. Our bound on AC does not, for that rea- 
son, depend on any conjecture concerning the size of EW 
boxes. Notice that this is not the case in the analysis 
of the "Magnificent Seven" [12, 81], which relies on the 
assumption AB'^'^ = AS'^' = to derive a probability 
distribution function for AC. 

The constraints on ACf^' and AC within CMFV fol- 
lowing from the simultaneous use of Rl Ab, A°^^, B{B ^ 
Xs'y), and B{B Xsl~^l~) can be seen in Fig. 4. All pan- 
els show frequentist 1 — CL levels. We see from the top 
and the lower left plots that two regions, resembling the 
two possible signs of Cf^{mi,), satisfy all existing experi- 
mental bounds. The best fit value for AC^** is very close 
to the SM point residing in the origin, while the wrong- 
sign solution located on the right in the upper (lower) 
panel is barely (not) accessible at 68% probabihty, as 
it c;orrcsponds to a B{B Xgl^l^) value considerably 
higher than the measurements [82]. In the upper (lower) 
panel of Fig. 4 the contribution Ai?'^' is scanned in the 
range ±0.1 (set to zero). In the former case the full re- 
sults read 



:>l+l- 



ACf = -0.039 ± 0.043 U 0.931 ± 0.016 (68% CL) , 
AC?*^ = [-0.104,0.026] U [0.874,0.988] (95% CL), 

while in the latter one we obtain 

ACf = -0.039 ± 0.043 (68% CL) , 



(10) 



ACf 



-0.104,0.026] U [0.890,0.968] (95% CL). 



(11) 



Similar bounds have been presented previously in [12]. A 
comparison of Eq. (10) with Eq. (11) makes clear that the 
size of AS'^' has only a moderate impact on the acces- 
sibility of the non-SM solution of ACy^ while it leaves 
the ranges themselves almost unchanged. Nevertheless, 
for \AB'*'' I > \B^s^' I the wrong-sign case AC^^' ~ 0.93 
cannot be excluded on the basis of B{B Xsj) and 
B{B Xsl'^l~) measurements alone. The same state- 
ments apply to AD although its impact on the obtained 
results is less pronounced than the one of AB' ' . No- 
tice that since the SM prediction of B{B Xs'f) [29] is 
now lower than the experimental world average by more 
than 1 a, extensions of the SM that predict a suppression 
of the b ^ sj amplitude are strongly constrained. In 
particular, even the SM point AC^''^ = is almost dis- 
favored at 68% CL by the global fit. This tension is not 
yet significant, but could become compelling once the ex- 



perimental and/or theoretical error on B{B — > Xg'j) has 
been further reduced. 

As can be seen from the top and the lower right plots 
in Fig. 4, in the case of AC only small deviations from the 
SM are compatible with the data. In the upper (lower) 
panel of Fig. 4 the contribution AS' ' is varied in the 
range ±0.1 (set to zero). In the former case we find the 
following bounds 



AC = -0.037 ± 0.266 (68% CL) 
AC = [-0.493, 0.387] (95% CL) 



while in the latter one we get 



AC = -0.026 ± 0.264 (68% CL) , 
AC = [-0.483, 0.368] (95% CL) . 



(12) 



(13) 



These results imply that large negative contributions to 
C that would reverse the sign of the SM Z-penguin am- 
plitude are highly disfavored in CMFV scenarios due to 
the strong constraint from the Z ^ bb POs, most no- 
tably, the one from We stress that we could not have 
come to this conclusion by considering only flavor con- 
straints, as done in [12], since at present a combination 
of B{B A57), B{B -> Xsl+l-), and B{K+ ^ tt+i^i?) 
does not allow one to distinguish the SM solution AC — 
from the wrong-sign case AC « —2. Eqs. (12) and (13) 
also show that the derived bound on AC is largely in- 
sensitive to the size of potential EW box contributions 
which is not the case if the Z ^ bb POs constraints are 
replaced by the one stemming from B{K^ — > -k^vu). 

It is easy to verify that the derived 95% CL bound 
for AC holds in each CMFV model discussed here. By 
explicit calculations we find that the allowed range for 
AC is [0,0.12] and [0,0.13] in the THDM type I and 
II [64], [-0.05,0.06] in the MSSM [65], [0,0.04] in the 
mlJED scenario [52], and [-0.07,0] in the LHT model 
[60] with degenerate mirror fermions. 

Other theoretical clean observables that are sensitive 
to the magnitude and sign of AC are the forward- 
backward and energy asymmetries in inclusive and exclu- 
sive b — > sZ+Z~ decays [79, 83] and the branching ratios 
of Kl TT^vi> and K+ tt+vD [12, 84]. 

BaBar and Belle have recently reported measurements 
of the forward-backward asymmetry Ay^{B — > K*l'^l~) 
[85, 86]. Both collaborations conclude that NP scenarios 
in which the relative sign of the product of the effective 
Wilson coefficients C|'^(mb) and CJq ("^b) is opposite to 
that of the SM are disfavored at 95% CL. While these 
results also point towards the exclusion of a large de- 
structive NP contribution to the Z-penguin amplitude, 
it is easy to verify that the present Ay^[B K*l~^l~) 
constraint is less restrictive^^ than the existing data on 



Assuming AS'^' = 0, AD < 4, and neglecting all theoretical 
uncertainties leads to the very crude estimate AC ~ [—1, 1]. 
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FIG. 5: The upper (lower) panel displays future constraints on 
ACf and AC within CMFV scanning AB'"' in the range ±0.4 
(set to zero) that are based on flavor observables only. The colors 
encode the frequentist 1 — CL level and the corresponding 68% 
and 95% probability regions as indicated by the bars on the right 
side of the panels. See text for details. 



Z bb considered by us. Notice that a combination of 
the branching ratios and forward-backward asymmetries 
of inclusive and exclusive b —^ sl'^l~ transitions might in 
principle also allow one to constrain the size of the NP 
contributions AS' ' and AD. A detailed study of the 



-> s/+Z measurements on 
AC, and AD is however 
beyond the scope of this article. 

The remarkable power of the POs 



impact of all the available b - 
the allowed range of Ai?'^' , 



Ah, and A 



Q.b 
FB 



in unraveling NP contributions to the universal Inami- 
Lim function C is probably best illustrated by a com- 
parison to one of the undisputed "heavyweight champi- 
ons" in this category, the K'^ — > n^vD branching ra- 
tio. A careful analysis of this decay shows that even 
under the hypothesis AB'^'^ = a future measurement of 
B{K^ TT^vv) close to its SM prediction with an accu- 
racy of ±10% would only lead to a slightly better 95% CL 
bound than the one given in Eq. (12), while already rel- 
atively small deviations of B'^'^ from its SM value would 
give the Z ^ bb POs the edge. This feature is illus- 
trated in Fig. 5 which shows the constraints on AC and 
AC^*^ following from the simultaneous use of the present 
B{B Xs^) and B(B Xsl^l~) world averages and a 
future accurate measurement of K'^ tt'^vv leading to 
B{K+ TT+vu) = (7.63 ±0.76) x 10"". In the upper 
(lower) panel the contribution AB^'^ is allowed to float 
freely in the range ±0.4 (set to zero). In the former case 
we find the following AC bounds 



AC = -0.057 ± 0.588 (68% CL), 
AC = [-0.768, 0.668] (95% CL) , 

while in the latter one we arrive at 

AC = 0.026 ± 0.282 (68% CL), 
AC = [-0.376, 0.420] (95% CL) . 



(14) 



(15) 



However, this "upset" of the tt^vv mode should 

not be over emphasized. While in MFV models the rates 
of the rare K — > ■kvu decays can be enhanced only mod- 
erately [12, 32] a very different picture can emerge in 
non-MFV scenarios with new sources of flavor and CP 
violation. Since now the hard GIM mechanism present 
in the MFV s — > dvv decay amplitude is in general no 
longer active, large departures from the SM predictions 
are still possible without violating any existing experi- 
mental constraint [60, 87]. Precise measurements of the 
processes —> ir'^vv and Kj^ — s- tt'^i/P will therefore 
have a non-trivial impact on our understanding of the 
flavor structure and CP violation of NP well above the 
TeV scale. This statement remains true even after tak- 
ing into account possible future constraints on the mass 
spectrum obtained at the LHC and the refinement of the 
flavor constraints expected from the B-factories [32, 60]. 



r 



In order to allow a better comparison with the re- 
sults presented previously in [12], we will set AB'^^ = 
when determining the allowed ranges for 



the branching ratios of K'^ tt+i^P, Kl n^vv, 
Ki^ fi~^ fi~ , B XiigVp, and Bj^s — > within 
CMFV. The corresponding lower and upper bounds at 
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TABLE III: Bounds for various rare decays in CMFV models at 95% probability, the corresponding values in the SM at 68% and 
95% CL, and the available experimental information. See text for details. 





Observable 


CMl V (957c CLj 


bM (68% CL) 


bM (95% CL) 


Experiment 






f4 29 10 721 


7 15 + 1 28 


15 40 9 111 


(14 7+i:^-<n [701 


B{Kl 


n°w) X 10" 


[1.55,4.38] 


2.79 ±0.31 


[2.21,3.45] 


< 2.1 X 10" (90% CL) [88] 


B{Kl ■ 


-> M+/U")SD X 10^ 


[0.30, 1.22] 


0.70 ±0.11 


[0.54, 0.88] 




B{B 


Xauu) X lO'' 


[0.77, 2.00] 


1.34 ±0.05 


[1.24, 1.45] 




B{B 


^ Xsw) X 10^ 


[1.88,4.86] 


3.27 ±0.11 


[3.06,3.48] 


< 64 (90% CL) [89] 


B{Ba 


X IQio 


[0.36, 2.03] 


1.06 ±0.16 


[0.87, 1.27] 


< 3.0 X 10^ (95% CL) [90] 


B{Bs 


H+fj,-) X 10^ 


[1.17,6.67] 


3.51 ±0.50 


[2.92,4.13] 


< 9.3 X 10^ (95% CL) [91] 



95% probability arc reported in Tab. III. For compari- 
son, we also show the 68% and 95% CL limits in the SM, 
obtained using the CKM parameters from a standard UT 
analysis. The calculations of the SM branching ratios all 
employ the results of [92]. In addition we take into ac- 
count the recent theoretical developments of [93-95] in 
the case of — !■ ir'^vv and Kl tt^vv, and of [96] for 
what concerns Kj^ A*"""/^"- contrast to the standard 
approach we normalize the B Xd^svv decay width to 
the B — > Xy^ev rate, while wc follow [97] in the case of 
Bd,a ~^ 1^^ 1^^ 1 since both procedures lead to a reduction 
of theoretical uncertainties. The actual numerical analy- 
sis is performed with a modified version of the CKMfitter 
code. The used input parameters are given in App. B. 

It is evident from Tab. Ill that the strong bound on 
AC, coming mainly from the existing precision measure- 
ments of the Z — > 66 POs, does only allow for CMFV de- 
partures relative to the SM branching ratios that range 
from around ±20% to at most ±60% for the given rare 
K- and B-decays. While our upper bounds are in good 
agreement with the results of [12], the derived lower 
bounds are one of the new results of our article. A strong 
violation of any of the 95% CL bounds on the consid- 
ered branching ratios by future measurements will im- 
ply a failure of the CMFV assumption, signaling either 
the presence of new effective operators and/or new fla- 
vor and CP violation. A way to evade the given limits 
is the presence of sizable corrections ^Cnp and/or AB""^ 
and AB' Whik- these possibilit ies cannot be fuUv 
excluded, general arguments and explicit calculations in- 
dicate that they are both difficult to realize in the CMFV 
framework. 



V. CONCLUSIONS 

To conclude, we have pointed out that large contri- 
butions to the universal Inami-Lim function C in con- 
strained minimal-flavor-violation that would reverse the 
sign of the standard Z-penguin amplitude are highly dis- 
favored by the existing measurements of the pseudo ob- 



servables Ab, and Ap^ performed at LEP and SLC. 
This underscores the outstanding role of electroweak pre- 
cision tests in guiding us toward the right theory and 
immediately raises the question: What else can flavor 
physics learn from the high-energy frontier? 
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APPENDIX A: VERTEX FUNCTIONS 

Below we present the analytic expressions for the non- 
universal contributions to the renormalized LH Zdjdi 
vertex functions in the CMFV models considered in this 
article. All expressions correspond to the limit of on- 
shell external quarks with vanishing mass and have been 
obtained in the 't Hooft-Feynman gauge. 

In the THDM type I, the only additional non-universal 
contribution to the Z — > djdi transitions stems from dia- 
grams with charged Higgs boson, , and top quark, t, 
exchange. An example of such a graph is shown on the 
top left-hand side of Fig. 2. The correction depends on 
the mass of the charged Higgs boson, M^, the top quark 
mass, m,i, and on the ratio of the vacuum expectation 
value of the Higgs doublets, tan/3. Using the decomposi- 
tion of Eq. (7) we find for the corresponding form factor 
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■THDM 



2- tan^ /3 



(224 - 9)A^g' - 3 (64 - ^Wt 



2 



+ 



96 Ml 

2 (2s^^ - 3) (2Mg^ - (4m^ + g^)M^^ + 2(mf + g^)m|) , 
^o(-Mjj j 



(M± 



±2 



- 6(24-l)(2My-2mra-)^^(^2^^.2^^.2) 



+ 



+ 



+ 



8s^^ (2Mg^ - 2m^ - g^) ^ , 2 2 2, 
^ '-Bo{q^,m;,mt) 



12(24-1) (M±4_2M±W 



to; 



2„2^ 



TO^q 



■Co(g^O,0,M±^M|^TO2) 



4 (45^^M^-^ + Aslmt - {8slM^^ - {Asj - 3)g^) 



-Co((Z^0,0,TO2,TO2,M±2) 



(Al) 



Here and in the following the coefficients AQ{m'^), 
Bq{p'^, ml, ml), and C'q{pI,pI, {P1+P2Y ,m\, ml, ml) de- 
note the finite parts of the scalar one-, two-, and three- 
point functions in the MS scheme as implemented in 
LoopTools [26] and FF [27]. The above result agrees with 
the findings of [28]. In the case of the THDM type II the 
overall factor tan^ /3 has to be replaced by cot^ /?. 



In the case of the MSSM only SUSY diagrams involving 

c;hargino. . aud up-squark, Ua, exchange lead to non- 
universal correction to the renormalized LH Zdjdi vertex 
of Eq. (7). An example of a possible contribution can be 
seen on the top right side of Fig. 2. The corresponding 
form factor can be written as 



C^±{q') = - 



AMI 



^ iX^-'')j,{X^-)a. \ 2M± M±,Co(?2, 0, 0, M±f , M±f ^ m|j U^iK^Sat 



m,n=l a, 6—1 



+ 



ln(M|;) + ^ ( + 2^o(M|;) + 2Ao(M|^^) - 4^o(m^ J 



M. 



±2 



2ml+q')Bo{q\M^^,M^^) 



-2^M 

+ 4 (M|^ - miJ{Mf^ - ml) Co(g^ 0, 0, M|;^ M|^, mi) 



- ln(M±f ) + -^[q'+ 4Ao(M?;) - 2Ao{ml) - 2Ao(m^ J 



2g2 



±2\ 



2M±f - ml 



m,~ 



q^j ^o(«^"^? ,TO? ) 



4 (M±f + mlml - {ml + ml - q') M±f ) Co{q', 0, 0, m?^, m?,, M±f ) 



+ 



-(224-9)M|;' + 3(64-1)to|^ 
6(M±J-m?J 

{2sl - 3) {2Ml* + 2to4^ - (4TOi - q^) M^^) 
3{Mf - ml )2g2 
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{2sl - 3) (2M±f + 2mt - 2 (2m|^ - q^) M±f - m?^g2^ 
+ 3(M--m|JV ^"^"^'^"^ 



+ 



44(2Mj;^-2m|^ + g- 
2 (2s2^ - 1) (M±f + - (2m|^ - q') M? 



2 



85^ (M±f + ml^ - {2ml ' 



3g2 



(A2) 



where = {8V2GFe^V;jVu)-^ with V being the CKM 
matrix and i.j = d,s,b. 

The LH charginoup-squark-down-quark coupling ma- 
trix takes the form 

x^^ = I w*ir^^ - w*2T^^ I y . (A3) 

Here and in the following sp = sin/3, cp = cos/3, = 
tan/3, etc. 

The unitary mixing matrices U and W are defined 
through 

C/*M^± = diag(M± , M± ) , (A4) 

with M^^ ^ being the physical chargino masses that sat- 
isfy M^^ < . M^± denotes the chargino mass matrix. 



which in terms of the wino, M2, and higgsino mass pa- 
rameter, /U, reads 

M,. = ( ^^M„s,\ 
The 6x3 matrices 

{T^-U = (r^)« , (r^«)„i = (r^)„,,+3 . (A6) 

are building blocks of the unitary matrix that diago- 
nalizes the 6x6 mass-squared matrix M?. of the up-type 
squarks: 

pc/t ^ diag(m?^ , • . • , m? J . (A7) 
In the super-CKM basis [99] , M? is given by 



u 



(A8) 



where Mfy^ = ttiq^ 1 and M^^ = rriuj^ 1 are the left and 
right soft SUSY breaking up-type squark mass matrices, 
Mu = diag(m„,mc,mt), Au = A^t contains the trilin- 
ear parameters and 1 represents the 3x3 unit matrix. 
We assume CP conservation, so all soft SUSY breaking 
terms are real. 

The result given in Eq. (A2) agrees with the one of [34]. 
To verify the consistency of the results one has to take 
into account that in the last equation of [34] the coeffi- 
cient Cii should read C12, and that arbitrary constant 
terms can be added to the second and third coefficient of 
the four different coupling structures in Eq. (A2) , since 
their contribution disappears after the summations over 



m, n, a, and b have been performed. In addition, the ex- 
plicit ln(M^^^) terms are absent in [34]. They have been 
chosen such that C^±{q^) coincides for = with the 
expression for the one-loop Z-penguin function given in 
[64]. 

In the mUED model diagrams containing infinite tow- 
ers of the KK modes corresponding to the M^-boson, 
W'^y the pseudo Goldstone boson, G^^, the SU{2) 
quark doublets, Qg(fc), and the SU{2) quark singlets, 
Uqi^k), as well as the charged scalar, a^^, contribute to 

the non- universal correction to the Zdjdi vertex. A pos- 
sible diagram can be seen on the lower left side in Fig. 2. 
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The only additional parameter entering the form factor pactification radius 1/R. We obtain 
in Eq. (7) relative to the SM is the inverse of the com- 



a 



mUED 



96 ^ 

k=l 



(-8(2^ 



3)m2 



(ft) 



{Usl - 27)Ml + 3 (6s2 
(m2 - Ml)Ml 



l)m( ) m: 



2(2^-3) 



+ 



2 (2s2„ - 3) ((2m| - (4M^ + 3^2)^2 + 2Mt + 7MW) Ml + (TM^ - 3m2) m^^^q^^ 



-^0(M^(fc)) 



2 {2sl - 3) (6M4 + im\W + SM^qA 
+ '-^^2 '-Mmfk,) 

2 {Sslmj + 2 ((85^ - 9)Mg, + 2sW) - {isj - 3)(6Mg, + 5q^)Ml) ^ ^ , ^ ^ 
MT^2 ,rn^^k),mt^k)) 

, 6 {{4sl - 2) m2 + 2(4.^ - 5)Ml + {2sl - 1) g^) ^2 ^ 

+ ^72^^2 ^0^1 ^^wikY^wik)) 

2(44-3)(6M2 +5g2) 

^2 ^0(a >TO(fc)>"^(fe)j 

+ + ((44 - ml + (44 - 3) q") mt 



- 3) (3M4 +3m^,)g2 ^4M^g2 ^ 2^4^ M^)Co(5',0,0,m,%),m,%),M^(,)) 



+ ((2«'- - 1)^* + (2 (4 - ml + (24 - 1) 9') 

- {^{Wsl - 11)A4 - (2^2^ - l)ml,y + 2 {si + l)Mlq') 

- 24 (3Mt + Sml^q^ + AMlq^) Ml) Co(g^ 0, 0, M^(,), M^(fc), m,%)) 

244 (3M^ + 3m2 ^ 4M^g2^ 

+ ^ ^2^^ Co(9^ 0, 0, M^(,) , M^(,) , m^,)) 

4 (4s^ - 3) (3M4 + 3m2 ^ + 2q') 

^ Co{q', 0, 0, m^,), m^,), iV4(,)) 



(A9) 



where rut^k) = ^J^t+n^fky ^w(k) = \J^l + mlky 
and = k/R. We note that our new result for 

CmUED((7^) coincides for _ q ^ji^jj ^^jjg one-loop KK 
contribution to the Z-penguin function found in [52]. 

In the case of the LHT with degenerate mirror fermions 
the only new particle that effects the Z djdi transi- 
tion in a non-universal way is a T-even heavy top, T+. 
A sample diagram involving such a fermion is given on 



the lower right-hand side of Fig. 2. The form factor de- 
pends on the mass of the heavy top, which is controlled 
by the size of the top Yukawa coupling, the symmetry 
breaking scale /, and the parameter xj^ = Xl/{Xl + X^)- 
Here Ai is the Yukawa coupling between t and T_(_ and 
A2 parametrizes the mass term of T+. Our result for the 
form factor entering Eq. (7) is given by 
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+ 



mh - ^t) (((6si - l)M|^ - (105^^ - 7)Ml) Ml - (64 - 1)(M|^ - M^)m|) 

{Ml^-Ml){ml-Ml)Ml 

2 {2sl - 3) (2M6^ + M4^g2 _ g(^2^ + q^)Ml^Ml + 4(M^ + 2(72)m4 ) 

2 (24 - 3) (2mf + - 6(M^ + q^) m|Mg, + 4(Mg, + 2q^)Mi) 
{ml - Ml)miq^ 

2 (24- - 3) (m2-M2) . 



-Ao(M^J 



Ao(m?) 



^2 ^^m^ 



+ 



- (4M^^ - i8Ml + q^)M^^ + 4(M^ + q^)Ml 
+ (2M|^ - 4(M^ + g2)M|^ + 2Mt + TM^g^) M^)Ao(M^) 
8s?, (2M4^ + {2Ml + g2)M2^ - 4M4 - 6M^g2 



Bo{q\M^^,M^^) 



8 (24m| + (2(5?, - ml + sW) ml - (24 - 3)(2Mg, + 3g2)Mg.) ^ 2, 



24 (mI^ + mf _ 2M^ - 3^^ j 



M2,g2 



Boiq^mlM^] 



6(M|^ - m2) ((4^2^ - 2)(M|^ + m?) + (44 - 6)M^ + (24 - l)^^) 



Bo{q\Ml,Ml) 



(AlO) 



164 (^T+ + M^y - (3M^ + 2q^)M}^Ml + 2{Ml + g2)2M^) 



MW 



Co(g^O,0,M|^,M|^,M^) 



+ lLf27^ (^^'^t - ml - (24 - 3)g2) - 2(3(s2^ - 2)Ml + 2(4 - 3)9^) m^M^ 
+ 2(2s^ - 3)(M2 + g2)2M^) Co(9', 0, 0, m?, m^, M^) 
■ ( {{2Ml + q')Ml 2{Ml + q^)Ml) m? 

2{Ml + q^) (m|^ -Ml- q^) Ml) Co(g^ 0, 0, ml , Ml) 
({2sl - 1)M|^ - {2Ml - {2sl - l)q^) M^^ - {{Qsl - 7)M^ + AsImW) 

- 4 4(M^ + 2q^)Ml) C'oiq^ 0, 0, Ml,Ml,M^^) 

+ ((2«- - 1) ™? - (2< - (2^'- - 1)9') ™* - ((64 - 7)M4 + 4s^M^«2) ^2 

- 44(M^ + 2g2)M4 ) CoiqlO, 0, M^, M^, m^) 



24 
MW 



12 

MW 



where Mt+ = f /vmt/ \/xl{1 — xl) and v ~ 246 GeV. work in 't Hooft-Feynman gauge is crucial for this com- 
Our new result for Clht(9^) resembles for g'^ = the parison. 
analytic expression of the one-loop correction to the low- 
energy Z-penguin function [60] . Taking into account that 
the latter result corresponds to unitary gauge while we 
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TABLE IV: Parameters that enter the standard and universal 
UT analysis. 



Parameter 


Value ± Error (s) 


Reference 


IKdl 


0.97377 ± 0.00027 


[30] 


\Vus\ 


0.2257 ± 0.0021 


[30] 


I^c6| 


(41.7 ±0.7) X IQ-^ 


[30] 


\Vub\ 


(4.31 ±0.30) X 10"^ 


[30] 


\eK\ 


(2.232 ± 0.007) X 10"^ 


[30] 


ArriK 


(3.4833 ± 0.0059) x 10"^^ MeV 


[30] 




(0.507 ± 0.004) ps^i 


[76] 




(17.77 ±0.10 ±0.07) ps"^ 


[100] 


sin(2/3)f,^ccs 


0.675 ± 0.026 


[76] 


mj(0 


(497.648 ± 0.022) MeV 


[30] 


ruBa 


(5.2793 ± 0.0007) GeV 


[30] 


rriBs 


(5.3696 ± 0.0024) GeV 


[30] 


£ 

JK 


(159.8 ± 1.5) MeV 


[30] 


-f d1/2 


(244 ± 26) McV 


ri r\-i 1 

[101] 




(281 ± 21) MoV 


[102] 


Bk 


0.79 ±0.04 ±0.09 


[103] 


fir. 


1 9s -1- n r\A 4- n no 
i.zo ^ u.u^ ^ u.uy 


[iU4J 


Bb, 


1.30 ±0.03 ±0.09 


[104] 




1 210+°°*^ 

^•^^'-'-0.035 


[101] 




0.5765 ± 0.0065 


[105, 106] 


net 


0.47 ± 0.04 


[106] 


Vcc 


1.56 ±0.37 


[106] 


Vb 


0.551 ± 0.007 


[105] 



APPENDIX B: NUMERICAL INPUTS 

In this appendix we collect the values of the experi- 
mental and theoretical parameters used in our numeri- 
cal analysis. The Higgs mass and the various renormal- 
ization scales are scanned independently in the ranges 
100 GeV < M° < 600 GeV, 100 GeV < Ht < 300 GeV, 
40 GeV < ^i„ < 160 GcV, 2.5 GeV < < 10 GeV, and 
1 GeV < /Uc < 3 GeV, respectively. 

The other parameters are displayed in Tabs. IV and V. 
Errors arc indicated only if varying a given parameter 
within its la range causes an effect larger than ±0.1% 
on the corresponding result. When two errors are given, 
the first is treated as a Gaussian 1 a error and the second 
as a theoretical uncertainty that is scanned in its range. 

Tab. IV contains the quantities that are relevant for 
the standard and universal UT analysis. We recall that 
there is a discrepancy of axound 1 a between the value of 
jV^bj obtained from inclusive and exclusive b — > ueP tran- 
sitions. Since \Vub\ only enters the standard UT analysis 
its actual value has no impact on our main results. We 
therefore use the weighted average of given in [30]. 
In the case of fs^ B^^ and ^ we take the values quoted in 
[101], which combines the values of the bag parameters 



TABLE V^ Quantities that are necessary for the calculation of 
the Z — > 66 POs and the rare and radiative K- and B-decays. 



Parameter 


Value ± Error(s) 


Reference 


Gf 


1.16637 X lO-'^ GeV-2 


[30] 


*W 


0.2324 ±0.0012 


[3] 




1/137.036 


[30] 




0.02768 ± 0.00022 


[107] 


as{Mz) 


0.1189 ±0.0020 


[30, 108] 


Mz 


(91.1875 ± 0.0021) GeV 


[3] 


Mw 


(80.405 ± 0.030) GeV 


[30] 


mt,pole 


(170.9 ± 1.8) GeV 


[109] 


ml' 


(4.68 ± 0.03) GeV 


[110] 


mc{mc) 


(1.224 ± 0.017 ± 0.054) GeV 


[111] 




105.66 MeV 


[30] 




0.1061 ±0.0017 


[112] 


C 


0.58 ±0.01 


[110] 


Ai 


(-0.27 ±0.04) GeV^ 


[110] 


A2 


0.12 GcV^ 


[30] 




(2.229 ±0.017) X 10""^ 


[95] 


K+ 


(5.168 ±0.025) X 10"" 


[95] 




(2.009 ± 0.017) X 10"^ 


[96] 


Aem 


-0.003 


[95] 


r{Bd) 


(1.527 ± 0.008) ps 


[76] 


r{Bs) 


(1.461 ± 0.040) ps 


[76] 



Bb^ determined by the JLQCD Collaboration using two 
light flavors of improved Wilson quarks [104] with the 
staggered three flavor results for the i?-meson decay con- 
stants obtained by the HPQCD Collaboration [113]. 
The central value and error of fs^Bg^ are taken from 
the recent publication [102] of the HPQCD Collabora- 
tion. For a critical discussion of hadronic uncertainties 
in the standard CKM fit we refer to [114]. 

Tab. V summarizes the remaining parameters that en- 
ter the determinations of the Z ^ bb POs and the rare 
and radiative K- and B-decays. In the case of as{Mz), 
we adopt the central value from [108], but rcscalc the cor- 
responding error by a factor of two to be consistent with 
[30]. We recall that the parameters kl, k;_|_, and scale 
like (A/0. 225)^ and that the values given in Tab. V corre- 
spond to A = \Vus\ — 0.225. This scaling has to be taken 
into account in order to find consistent results in the case 
of the rare JsT-decays. The IR finite long-distance QED 
correction factor Aem entering the prediction of 
TT^i/D accounts for photon emission with energies of up 
to 20 MeV [95]. Since mc{mc) and the phase-space factor 
C = |Kb/K6pr(B ^ X^ei>)/r{B ^ X,ei?) are strongly 
correlated we take both of their values from global analy- 
ses of semileptonic i?-decay spectra [110, 111]. However, 
to be conservative, we treat the errors of mdmc) and C 
as independent in our fit. If the anti-correlation would be 
included, the individual uncertainties from m,c{mc) and 
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C in the branching ratios oi B ^ Xs^, B Xgl^l , each other [29]. 
and B X^^gvv would cancel to a large extent against 
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